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Abstract. We consider non-local currents in the context of quantized afSne algebras, following the con- 
struction introduced by Bernard and Feldcr. In the case of Uq{A^^-') and Uq{A^^), these currents can be 
identified with configurations in the six-vertex and Izergin— Korepin nineteen-vertex models. Mapping these 
to their corresponding Temper ley-Lieb loop models, we directly identify non-local currents with discretely 
holomorphic loop observables. In particular, we show that the bulk discrete holomorphicity relation and its 
recently derived boundary analogue are equivalent to conservation laws for non-local currents. 



1. Introduction 

Discretely holomorphic observables are correlations functions in a two-dimensional lattice model which 
satisfy a discrete version of the Cauchy-Riemann (CR) equations. In the context of the Ising model, lattice 
fermions with this type of property were first described in [13) . More recently, the discrete CR equations were 
used by Smirnov as a basic tool to study rigorously the scaling properties of Ising interfaces |17] . They were 
then exploited in the Probability literature, to obtain mathematical proofs of several Coulomb-gas results. 
For instance, still in the Ising model, this approach yielded the scaling limit of domain walls and Fortuin- 
Kasteleyn cluster boundaries [101 US]: and the spin and energy correlation functions [9]. For self-avoiding 
walks, it provided a rigorous way to determine the bulk |14j and boundary [4] connectivity constants, and it 
has also proved very useful for numerical purposes [S]. 

In the meantime, some discretely holomorphic observables have been identified in other 2D lattice models, 
including the Z^r clock model [24^ and the dense ^5] and dilute [18] Temperley-Lieb (TL) loop models. 
These observables are essentially non-local, either because they include disorder operators (in spin models), 
or because they are defined in terms of an open path attached to the boundary (in loop models). In all 
the known examples, it was observed that the discrete holomorphicity condition is satisfied precisely when 
the Boltzmann weights are such that the model is integrable [211 HI]. Recently, this statement was also 
extended to the boundary Boltzmann weights [171 HI]- These observed relations between the notions of 
discrete holomorphicity and integrability have been explored further recently [I], but they still call for a 
more systematic understanding: this is the object of the present work. 

An obvious starting point for this is to try and construct discrete holomorphic observables from the 
underlying symmetries of a the lattice model. This idea is very reminiscent of the construction of non-local 
conserved currents ip{z) in lattice models possessing a quantum group symmetry [5]. Indeed, in that context, 
ipiz) is non-local because it includes a path connecting z to a reference point, in a similar way to disorder 
operators in spin systems. Moreover, the conservation equation for the current is a linear relation between 
the values of at the points adjacent to a given vertex, like the discrete holomorphicity condition. These 
resemblances jl5j between discrete holomorphic observables and conserved currents can actually be made 
more precise. 

The present paper is based on a simple observation: in the case of an affine quantum group symmetry, the 
conservation equation of currents can be written as a discrete holomorphicity condition on the rhombic lattice 
with opening angle a, provided an appropriate relation between a and the spectral parameter is introduced. 
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We thus consider two simple loop models where discretely holomorphic observables - which we shall call 
for short loop observables - are known: the dense and dilute Temperley-Lieb models on the square lattice. 
Using the mapping of these loop models onto vertex models possessing Uq{A^i ^) and UqiA':^ ^) symmetry 
respectively, we construct the conserved currents, and show that they map to the loop observables identified 
previously. This analysis is extended to boundary observables in the case of general diagonal integrable 
boundary conditions. 

This point of view explains the somewhat mysterious observation of |24[ 118) that discrete holomorphic- 
ity somehow "linearizes" the Yang-Baxter equation by providing us with a linear equation for integrable 
Boltzmann weights. Indeed, from our point of view, this linearization procedure is nothing but Jinibo's 
interpretation [20 of the i?-matrix of an integrable model as a representation of the universal i?-matrix of a 
quantized affine algebra, which by definition satisfies such a linear relation, as will be explained below. 

The plan of the paper is as follows. In §[21 we review the construction of conserved currents introduced in [6] 
and expose a general identity for the adjoint action in this context. §[3] reviews the correspondence between 
the six- and nineteen- vertex models and the dense and dilute Temperley-Lieb loop models respectively. In 
particular, we show that the integrable weights can be obtained by solving the intertwining relations, which 
are linear equations in the Boltzmann weights assigned to loop model tiles. In § 21 we use the mapping 
between vertex and loop models to express the currents as loop observables, satisfying the discrete Cauchy- 
Riemann equations. §[5]extends the work of [6] to systems with a boundary, and focuses on the interpretation 
of current conservation at the boundary. § [HI introduces boundary tiles into the dense and dilute Temperley- 
Lieb loop models. Integrable weights for these tiles are obtained by solving linear equations which are the 
boundary analogue of the intertwining relations. In § [71 we express our currents (which satisfy current 
conservation at the boundary) as loop observables obeying boundary discrete Cauchy-Riemann equations. 
In § [HI we use the Coulomb-gas approach to present the continuum limit of the loop observables. We conclude 
in §[3 



2. Vertex models, currents and quantized affine algebras 

2.1. Vertex models. In this section, we recall how vertex models in statistical mechanics can be defined in 
terms of Boltzmann weights which are given as intertwiners of representations of quantized affine algebras. 
More specifically, we will consider the case in which we have a quantized afltine algebra U with a spectral pa- 
rameter dependent module 14. We denote the associated representation as (tt^, 14) where tTz : U ^ End(T4) 
and additionally use the notation nz-^^...^ZM — ""zi "X) • ■ • (S) tt^^ for tensor products of such representations. 

Assuming that 14 ® 14) is generically irreducible, the i?-matrix R{z/w) : 14 fX) 14i — > 14 "Xi Ko that defines 
the vertex model is given as the solution (unique up to overall normalization) of the linear relation 

(1) R{z/w) 7r,,^(A(X)) = ^,,„(A'(X)) R{z/w) 

for all X ^ U; where A : [/ ^> [/ C/ is the coproduct, A(X) = ^ Xi X2, and A' the coproduct with the 
order of the tensor product reversed, A'(X) = ^ X2 <S^ Xi. We represent the i?-matrix pictorially by 

R{z/w) = z — > 

w 

The arrows drawn on the lines are purely to indicate "time fiow" : reading an equation from right to left 
corresponds to reading along a line in the direction of the arrows. 
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Let us define the multiple coproduct A^^) : U U®^^^ by A^^+i) = (A (g) l)A(^) and A^^) = A. The 
monodromy matrix {z;wi,..., wl) ■■ Vz <E) Vnji (E) ■ ■ ■ (E) K-i, -> 14 «) K-^ (g) ■ • • ® is defined as 

T^^\z;wi, . . .,wl) = Rql{z/wl) . . .i?oi(z/wi) , 

where the subscripts on i?- matrices indicate the evaluation modules in which they act, i.e., O Vz, 
1, . . . , i ^ 7 • • • , ■ Its graphical representation is 

T^^\z;wi,...,wl) = z 



Wl W2 ■ ■ ■ Wl 

Until specified otherwise, all lines will be oriented up/right in what follows, so that we shall omit arrows on 
pictures. 

By vertically concatenating M monodromy matrices T^^\zi;wi, . . . ,wl), 1 < * < M, one obtains a 
rectangular lattice il of width L and height M. Each horizontal (resp. vertical) line of this lattice is oriented 
from left to right (resp. bottom to top), and denotes the vector space Vz^ (resp. V^j). 

To simplify the discussion, from now on we will work on a horizontally and vertically homogeneous 
lattice, by assuming that all horizontal (resp. vertical) lines carry the same evaluation parameter Zi = z 
(resp. Wj = w). We will be interested in operators that act on the vector spaces encoded by the lattice, 
which graphically corresponds to the insertion of a "node" at an edge. The edges of the lattice will be 
denoted by coordinates pairs {x ± ■^,t) or {x,t ± i), where {x,t) is a vertex of 51, and x (resp. t) is the 
horizontal (resp. vertical) coordinate. 

Typically, one is interested in the case where the left/right boundaries of this lattice are fixed in some 
way. In such a case, each row of the lattice becomes an operator in V^i ® ■ • • g) Ko^ :— Vi (E) ■ ■ ■ ® Vl that 
we will rather loosely call the "one-row transfer matrix" T. We will use this construction below, despite the 
fact that for now we do not discuss the boundary conditions explicitly. 



2.2. Hopf algebras and graphical relations. Following Bernard and Felder [6], we consider a set of 
elements {Jo, Qj', Q°'b}, a,b — 1, . . . ,n, oi a. Hopf algebra U. The elements Qj' and G^fc are assumed to be 
inverses of each other: 

(2) e^t'e-, = 5a,c and e^Gfc" = <5a,c 

(where here and subsequently repeated indices are summed over) while the coproduct A and antipode S of 
all elements have the form: 

(3a) A( J,) = J, 1 + Ga' ® Jb SiJa) - -e^ Jb 

(3b) A(ea^) = Qa' E> Ob' SiQa") = a 

(3c) A{e\)^e\(g>e\ 5(e°6) = efa". 

It is also useful to define Ja :— —Q^'aJb, which has the coproduct and antipode 

(3d) A{Ja) ^ JbE)e\ + l^Ja SiJa) = Q^JcQa''. 

Given two elements in { Ja}, say Ji and J2, Ji + J2 can be trivially added to the set { Ja} in such a way that 
(PatiSdp still hold. It is a little less obvious that the same is true of Ji J2, with appropriately defined sets 
{0a^} and {0°(,}- Therefore we only need to specify a set { Ja} that generates J7 as a unital algebra. 
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The i?-matrix, R : U ®U U (E) U, switches the order of tensor products in the coproduct: namely, 
RA{Xa) — A'{Xa)R for all Xa E U. Applying this to the coproducts in (|3aH3dl) gives, respectively 



(4a) 
(4b) 
(4c) 
(4d) 



R{Ja ® 1 + Oa^ Jfc) = (1 ® Ja + Jb «) ea'')R 

R{ea''<E>Qbn = (eb^®e/)i? 
i?(e"fc®e''c) = i&c®e\)R 

R{Jb «) + 1 ® Ja) - ®Jb + Ja® l)R . 



Suppose now that we have a representation (tt, V) of the Hopf algebra U. In the spirit of [6], we can represent 
7r(Ja), Tr{Qa^) and 7r(0°;,) by the following pictures (from now on we always discuss representations of U, 
and so suppress the appearance of the tt) : 



Ja= -^1, = 




where the vertical line denotes the vector space V with an upward arrow that we have suppressed, and 
subscripts (resp. superscripts) correspond to incoming (resp. outgoing) arrows. The operator Ja has the 
graphical representation 



Ja = — a that we simplify to Ja = B^-^ a 



so that a connection of a wavy line to a solid line from the right (compared to the direction of time) 
corresponds to the insertion of the operator Ja- 

Using these notations, the equations listed have natural graphical meanings. Relation © is expressed 
graphically by 



and the first coproduct relation (j4a|) is equivalent to 



(5a) 



a < 










1 



a 


a 








h 





R{Ja ® 1) 



RiQa" ® Jb) 



(1 ® Ja)R 



{Jb (E> ea'')R 
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where we recall that the "time flow" is south-west to north-east. Similarly, for "tail operators" one has 




which means one can move the tail freely across vertices. The remaining equations (|4cH4d| have analogous 
the graphical equivalents, but the tails now enter from the right: 





















a 


a 




2.3. Non-local currents and conservation laws in the bulk. Continuing along the lines of we 
act with the repeated coproduct A^^^ on the elements of U to define non-local currents. By iterating the 
coproduct in (|3a|. we obtain 

L 
x = l 

The object thus constructed, Jo, is the charge associated with the time component ja\x) of a non-local 
current. Acting on a tensor product Vi Vl, we have the graphical representation 




X 

where each solid line corresponds to a space Vi (the tensor product is ordered from left to right, as in 
the corresponding algebraic expression). If 1^ ~ as below, then each solid line will carry an index in 
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{1, . . . , d}, and each wavy line an index a G {1, . . . , n}. The intersection of a wavy hne and a sohd hne is a 
0a; acting on the sohd hne. 

We want to reinterpret (|4a| (or its graphical equivalent, ([5a|) as a discrete current conservation for a 
vector field {j'a\ja^), where the index x corresponds to horizontal direction and the index t to vertical 
direction. This leads us to define j^\x), x G Z + 1/2, as the insertion of a dot on the horizontal edge 
[x — 1/2, a; + 1/2] (with a tail attached, one half-step up and then to the left). If one wants to define jl"^^ as 
an operator ouVi® ■ ■ ■ ®Vl, one needs to "embed" it inside a transfer matrix: i.e., letting T be the one-row 
transfer matrix (recall that the boundary conditions on the left/right are left undetermined in this section), 
we define 

T^'^jI^\x)T^'^ = 



a 
























X 









Adding a tail to all terms in ([5al) which extends all the way to the left, we can straighten it using (|5b|) . 
Assuming that the tail commutes with the left boundary, we find 

(7) j^) (a; - 1/2, t) - j^) (x -I- 1/2, t) + J (*) (x, t - 1/2) - (x, t + 1/2) = 

where in the operator formalism, the time evolution of any operator O is given by 0{t) = T*C'T~*. Equation 
([7|) expresses the conservation of the current {ji^\ja'') that we have defined. 

Summing (jT)) over x results in the conservation law for the associated global charge 3a (l6|) , up to boundary 
terms: 

(8) 3a{t + 1/2) -3a{t- 1/2) = (1/2, t) - j^) {L + 1/2, t) . 
If we depict the charge 3a by 



where the oriented dashed line represents summation over the position of the node {i.e., discrete integration), 
then ([8]) has the graphical interpretation 





















> 








1 




-->-- 







In this paper, we shall be mostly concerned with the "local" relation ([7]) and not so much with the global 
relation (|5]). Note however that even the former relation is not strictly local because of the tails; we therefore 
shall have to pay attention to the boundary conditions to the left when applying it in what follows. 

It is the purpose of this paper to relate ([7]) to the so-called discrete holomorphicity condition. 



2.4. Adjoint action. A detailed discussion of the action of U on local fields and of its adjoint action can 
be found in §2.4 and 2.5 of [j^. Here we summarize some relevant facts. Let us write the general coproduct 
as 

A(X)=^X(1)®X(2), 
for any X G U. The adjoint action of a Hopf algebra is defined by 



adx (Y) :=^X(i) r 
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In the case of elements of the form Ja , whose coproduct and antipode are given by pap , this means that 

adj„ ( Jfa) = JaJb - Qa'JbQ'^cJd ■ 

The natural action of J a on Jfc, viewed as an operator on Vi (g) • • • Vl, is obtained by applying A*^-^); 



where Aa on the r.h.s. is the action of Ja on local fields, and is best described graphically: 



(9) 



= b 




Once again, dashed lines denote discrete contour integration. 



(x) 

A similar procedure can be carried out for the action on the other component of currents. In this 
way, one can organize various currents as multiplets of U (or of subalgebras oi U). In general, starting from 
some generators Ja, the adjoint action of the whole of U produces a large subspace inside U. Note however 
that in the case of quantized affine algebras, to be discussed now, the Serre relations imply that the module 
generated by the adjoint action of a Uq{Ai) := Uq{sl2) subalgebra on some other Chevalley generator is 
finite-dimensional. 

We now apply the formalism of previous sections to particular quantized affine algebras. 



2.5. Quantized affine algebras. In this paper we always take ?7 to be a quantized affine algebra Uq{Q) 
corresponding to a rank 1 affine Lie algebra g. Let Aij denote the entries of the generalized Cartan matrix 
for Uq{Q), and let di be integers such that diAij = djAji, whose greatest common divisor is 1. Then the 
Chevalley presentation of is given in terms of the generators {Ei, Fi,Ti}, i e {0, 1}, satisfying the list 

of relationqj 



(10) 
(11) 

(12) 
(13) 



T,T- 



I- Aij 

EH' 
EH' 

fe=0 



where we use the notation 



1- A 
k 

l-A^ 



m 
n 



iE,y--^'^-^E,iE,)'^ = 0, 



(F,)i--^-~'=F,(F,)" = 0, 



q—m^ ^^m — n+1 _ ^— m+n— 1^ 

(q"-g"")...(g-g-i) 



The coproduct of these generators is taken to be 

A{E,) ^E,®l + T,®E,, A{F,) ^ F, Tr^ + 1 ® 



For a review of quantized affine algebras see \8\. 
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It is also useful for our purposes to introduce the modified generators Ei :— q'^'^TiFi, whose coproduct takes 
the more convenient form: 

A{Ei) =E,®l + T,®E,. 

The correspondence between these generators and the Hopf algebra elements J a, Qa^, Q°'b introduced in i) 12.21 
is immediate. For a = 0, 1 we set: Ja — Ea,Ea, their coproduct being of the form of (I5al) : Qj' — Sa.bTa 
(resp. 0";, ~ Sa.bT^^^) their coproduct being of the form (j3b|) (resp. ([5c|) ): Ja = Fa, their coproduct being of 
the form pd)l ). 

In the following subsections we describe the two quantized affine algebras which will interest us in this 
paper, as well as the representations to be considered. 



2.5.1. The case U = Uq{A\^^). The first case of interest to us is Uq{A[^^), for which the generalized Cartan 
matrix is given by 



Aio An 



2 -2 
-2 2 



and do = di = I. The representation (tt^, Vz) is taken as the level-zero fundamental (principal) evaluation 
representation Vz = C^[[2]] 



The i?-matrix R{z) is given by 



TTziEo) 

nz{Ei) 



z-i 







q-^ 
q 

q 
q-^ 



(14) 



Riz) 



\ 



qz — q ^z ^ 













z — z^^ 


q - q-^ 








q - q-^ 


z ~ z^^ 











qz - 





which gives the weights of the 6-vertex model in the principal gradation. 

2.5.2. The case U^UqiA'^ '). The second case which we study is Uq{A2 ), for which the generalized Cartan 
matrix is 



-^00 Aoi 
Aio An 



2 -4 
-1 2 



and do ~ 1, di ~ 4. The representation (tt^, Vz) is now Vz = C'^[[z, z^]] with 

T^ziEo 





z^-'^ifiiq) I 1 
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where ip{q) — {q + Q^^Y^^ , and £ is an arbitrary constant which we will fix later. The i?-matrix is 









Wl6 


CO5 


W8 




Wl9 


(15) 


R{z) = 






OJ12 


Wi 


Wi3 


We 








Wl8 




Wg 


a;4 


Wl7 


Wii 



V 



where the entries are given by 



WlO = Wii 



W2 = 

U)3 = U)5 

U)6 = Us 

UJ7 = iO(j 

^12 = 
Wl4 = Wis 
W16 = W17 
W18 
Wl9 



{z- z ^){q^z 

„2 „-2 



Wl5 / 



Z+\q^-q-^){q'z + q-'z-') 



+^z+'{q^ 

- 1 \ / „3 , 



g-^)(z-z-^) 



(z-z-^)((7'^z + g-'^z-^) 

„-2„-l\/„3., I „-3_-l 



(z - z~^){qz 



^z-')(ci'z + q-''z-') 



■q-'z-') 



= Z+^' 



{q^-q-%{q^+q-^)qz^-{q-q-')q- 
{q^ -q-^W + q-^)q-'z-^ + {q 



q-')q'] 



These coincide with the Boltzmann weights of the Izergin-Korepin 19-vertex model, up to factors of z^^ 
which come from the fact that we have not yet fixed the gradation. 

3. From vertex models to loop models 

In this section we review the loop/vertex model connection for Uq{A^^^) and {/^(Aj^^) models in the bulk. 
Until this stage, our analysis has been on a lattice with an arbitrary angle between horizontal and vertical 
lines. We now specify our domain further, by considering a rhombic lattice of definite angle a, which we will 
ultimately relate to the ratio of the spectral parameters z/w. We also draw the dual lattice using dotted 
lines: 



(16) 



R 




where the edges of an elementary "plaquette" (elementary rhombus of the dual lattice) are of unit length. 



3.1. The dense Temperley Lieb model and the Uq{A'"p) vertex model. The dense Tcmpcrlcy Lieb 
model is defined by assigning weights a and h to the following two local configurations of a rhombus with 
top-left lattice angle a: 
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a b 

and weight t to closed loops in a given lattice configuration C. Thus, C is assigned the Boltzmann weight 

(17) M/(C) = a^»(C') ^iv,(c)^ 

where Na{C) (resp. Nb{C)) is the number of plaquettes of weight a (resp. b), and Ni{C) is the number of 
closed loops in C. 

Let us now recall how 

(18) T^~^{q + q-'), q:=-e- 

We then associate the above configurations with operators A, B: ViSV^V^V (acting SW-NE), with 
V — C It) ® C 14,). This is done by dressing each configuration with arrows and reading off the associated 



these weights are related to those of the six- vertex model. Firstly, we identify 

1\ „ 2Tziu 



weight e^^^ associated with total turning angle 50 In this way we obtain 



(19) 



/ 1 

e- 




V 









^ 
























1 ) 



B 



/ 







g-2iiy(7r- 





V 












/ 



in the basis {|tt) ^ Iti) 1 14-t) : lii)}- The ^-matrix is then the linear combination of these operators, dressed 
by their respective weights: 

fa \ 



(20) 



R:^ PR^a A + b B ^ 













V 







2"^"a + e^"'^''-°'^b 

a J 



3.2. Six- vertex model weights from intertwining relations. It is well known that the commutant of the 
quantized algebra Uq{Ai) acting on a tensor product of two-dimensional representations is the Temperley- 
Lieb algebra [19]. This suggests that the two terms A and B in ([20| are intertwiners for the subalgebra 
(^Ei, Ei,Ti^ = Uq{Ai). Indeed, one can check that 



(21) 




(A(X))=7r^,,(A(X))^ 
(A(X)) =^^,,(A(X))B, 



X = El, Ei,Ti , 



provided the following relation holds between angle and spectral parameters: z/w = e 



Notice that the 



equation for X = Ti is automatically satisfied, since both A and B preserve the total magnetization. 

Moreover, imposing that R commute in a similar way with the action of Eq, Eq fixes the ratio a/b, say 

(22) a = qx — q^^x~^, b = x — , x := z/w. 

Substituting these values for the weights into (|20p . we recover the 6- vertex i?-matrix (|14p . 

Remark: The explicit connection of the operators (|19p with the usual generators of the Temperley-Lieb 
algebra is as follows. On a strip of width L, the Temperley-Lieb algebra with generators {gi, . . . , ^l} satisfies 
the list of relations 

(23) 



9j9j±i9j 
9]9k 



9j 
T9j 
9k9]i 



if \j - fc| > 1. 



We assume, as shown on the pictures, that the loop hnes enter/leave orthogonally to the sides of the rhombus. 
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The generator gj acts non-trivially on spaces at positions j and j + 1, and the weights for a plaquette at this 
position are encoded in the Tt-matrix Ttj.j+i = a 1 + b gj. Both the identity 1 and TL generator gj have well 
known graphical interpretations, as a 45 degree rotation and deformation of the tiles above into squares. 

Consequently, we expect that R and TZ are related by a simple gauge transformation (corresponding to 
the passage from principal gradation to homogeneous gradation), which is indeed the case; we find that 



TZ^U-^ RU' , 



since under this transformation the two terms in R become 

/ \ 



U-^ AU' U-^ B U' 



-q-^ 1 
1 -q {) 

V y 



where the second term is the well-known spin-i representation of gj. 



3.3. The dilute Temperley Lieb model and the Uq{A2 ) vertex model. The dilute Temperley-Lieb 
model (or 0{n) model) is defined by the plaquette configurations: 



J 



Ml 



U2 



with corresponding weights shown beneath the configurations. The Boltzmann weight of a configuration C 
is given by 

w{c) = t^'(^) ^f-('^^Kf-(^'^f-(''^^f-(^V^^(^) , 

where Na is the number of plaquettes of type a, and Ni is the number of closed loops. In a similar way to 
the dense case, we introduce the parameters 



q:=e 



Now we identify the configurations with operators T, [/{, Ui' , U2,U2 ,V' ,V" ,Wx,W2 -.V^V ^V®V, where 
V — = C It) © C |0) ® C 11). We dress lines in the plaquette with arrows and associate them with |t) or 
4,), we identify missing lines with |0), and collect an associated weight e'^"^ for the total turning angle 6. To 
save space, we will not write down the explicit form of the resulting operators. 
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The i?-matrix is, as before, the hnear combination of ah operators dressed by their Boltzmann weights. 
In the basis {|tt) , It 0) , lU) , |0 t) , |00) , |0 ;) , |;t) , 11 0) , |U)}, it is given by 



(24) R^tT + ui {U[ + [/{') + U2 {U'2 + t/2') +v{V' + V") +^11^1+ W2 W2 = 





































V 























{wi + W2e-^'"''')e^""^ 










W2 











V 

































t 

































V 











W2 
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3.4. Nineteen- vertex model weights from intertwining relations. In analogy with the dense TL case, 
we wish to identify the plaquette operators as intertwiners of the subalgebra {Ei, Ei,Ti) . We find that 

Y ^,,„(A(X)) = 7r^,,(A(X)) Y, X = £^1, ^1, Ti , F = T, {/{, U[' , U'^, U'^, V , V" , VFi, VF2 , 

provided we set {z/w)'^^ — e~^"^". Furthermore, imposing that 

i? 7r,,„(A(X)) = 7r^,,(A(X)) R X = Eo,Eo,To, 

determines the Boltzmann weights up to normalization: 

t^{x~ x'^^){q^x + q'^x^'^) + {q^ - q~'^){q^ + q~^) 
Ml = (q^ - q^^)iq^x + q^^x~^) 

U2 = i{q^ - q^'^){x - x^^) 

V = {x - x^^)iq^x + q^^x^'^) 

Wl = (q^x - q^'^x^^)(q^x + q^^x^^) 

W2 = [x - x^^){qx + q^^x^^) , 

where we have again set x z/w. Inserting these values for the weights into (|24p . we recover the 19-vertex 
i?- matrix (|15p . 

4. Discrete holomorphicity in the bulk 

The goal of the present section is to connect ([7]) to discrete holomorphicity. To do so, we must finally specify 
our boundary conditions. In all cases to be considered, we shall choose "reflecting" boundary conditions, in 
which adjacent edges on the boundary are paired by a common loop (or are empty, which is allowed in the 
dilute TL model). The only exception to this rule will be two (or one, in the dilute case) boundary edges 
from which an open, unpaired loop propagates. 

The reason for making this choice is that these are simple boundary conditions for which the boundary 
trivially commutes with the tail operators Tq,Ti, allowing us to apply Indeed, our results extend to any 
choice of boundaries for which the tail operators satisfy this property. In particular, we wish to emphasise 
that our results do not require integrable boundary conditions. 

The discussion of integrable boundaries, and the associated boundary discrete holomorphicity, is deferred 

to gSHZl 



4.1. Application to the dense loop model. 
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4.1.1. Loop observables associated to Eq and Eq. Let us now consider the insertion on the lattice of the 
current eo associated to the operator Eq. We want to translate this insertion into the language of loops. We 
consider a model in which all loops are closed except one open path 7 that connects two fixed boundary 
defects. In the vertex model, this is simply done by setting free boundary conditions for the two boundary 
defects. An example of such a configuration is shown below. 




The important observation is that if the marked edge on which Eq sits belongs to a closed loop, then the 
contribution is necessarily zero because of conservation of m(1) charge. Graphically, 

Eq = 0, and -^^>-^^^^o^ = . 

To be precise, we have included the tail in these equations, but in fact the tail has an irrelevant contribution, 
since it is diagonal in the evaluation representation we are using. In subsequent pictures, the node marking 
will always correspond to insertion of Eq, and the wavy line to Tq. We conclude that in order to have a 
non-zero contribution, the marked edge must lie on the open path 7, and that the latter, whose orientation 
we had left unspecified so far, must be incoming at both boundaries: 




additional factor e'^'^'^^i'^) where 9{C) is the angle spanned by the portion of 7 between a boundary entry 
point and the insertion of Eq (the red or green line) - the two angles are equal. In fact it is easy to see that 
6{C) = Trk{C) where k{C) is an integer whose parity is fixed by the relative locations of marked edge and 
boundary entry points of the arc (on the example fc(C) = 2). Finally, there is a contribution from the tail 
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of q^"^ (wavy line). It is easy to check that if k is even, red and green hnes each cross (algebraically) the 
wavy line k/2 times, whereas if k is odd the red line crosses it [fc/2j times and the green line \k/2] times. 
In both cases we find that this produces a factor q''^^\ Using q — 6*^(21/- 1) putting everything together 
we find 

(4*)(x,0)=| E W{C) e<*^-'^'^^\ 
where we use the notation (O) for the expectation value of an operator O, and Z is the partition function, 

(x) 

Now let us apply the same reasoning to Cg , for which a typical configuration is shown below: 



First we have a factor z. Then we have a factor from the angle e^*"^ where this time 9 = —a + kn where 
k has fixed parity (on the example, k = —1). Finally, the tail crosses algebraically the red line [fc/2j times, 
and the green line \k/2] times. In total we get e2*'^eg(»+")/'^ = e*(4i^-i)egi(2>.-i)a| q^. 

Note that z e^*''" — w. Therefore, if we define a function on the edges of the lattice 

_,je^J\x,t), (x,t) eZx (Z+i) 
^ ' \e»"e^"^(a;,0, (x, t) £ (Z + i) x Z , 

then we have 

c|(x,t)e7 

where 9{C) is the angle spanned by the open path 7 from one boundary to {x,t). 
Furthermore, the current conservation equation ([7]) for ja = cq is 

(25) e^''^ (x - 1/2, t) - el^^ {x + l/2,t) + e^*^ (x, i - 1 /2) - e[,*' (x, t + 1/2) = . 
Rewriting (j25l) in terms of 4io{x,t) we find 

(26) Mx, t - 1/2) + e*("-"Vo(x + 1/2, t) - e'('^-"Vo(x - 1/2, t) - 4>o{x, t + 1/2) - , 



■^If wc instead defined q = e^^^'^^'^^'l we would obtain the total factor e«{4f-i)Se*(^''+^)'*, leading to an observable which is 
discretely anti-holomorphic. There is no preferred definition for q, since exactly half of the observables are anti-holomorphic 
regardless of which choice we make. 
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which is precisely the discrete holomorphicity condition around a plaquette of the type: 

(x,t+l/2) 

{x-i/2,t) • • {x+i/2,t) {x,t) e . 

■ • 

{x,t-l/2) 

This equahty is vahd at the operator level, i.e., when inserted in an arbitrary correlation function (• • •). 

Note that 0o is exactly the lattice holomorphic observable identified in [25j [18]: we have shown here 
that this observable can be obtained by the construction of conserved currents associated to the Uq{A\^-') 
symmetry of [B]. 

Everything can be repeated with Eq instead of Eq. This leads to the conjugate observable 



Z 

c\{x,t)ej 

which therefore satisfies an antiholomorphicity condition. 

4.1.2. Loop observables associated to Ei and Ei. There is a simpler observable, obtained by considering Ei. 
We use the same type of boundary conditions as above; the whole discussion goes through, except for the 
following modifications. Compared to E'o, the arrows on the open path 7 are inverted, but the tail is inverted 
as well (it is made of q^'^')- The result is that the tail and the angle factor almost compensate; for e^*^ we 
find using the exact same reasoning that 

C|(x,t)G7 

Note that e~*^ — (—1)'" is independent of the configuration and simply measures the flux of 7 through 

the marked edge (i.e., the only configurations that contribute to (e^*'') are those for which the average 

(2;) 

flux is nonzero, that is, where the open path goes through the marked edge). Similarly, for e\ one finds 

^-2iv0^{0+a)/iT ^ ^-ie^i{2v-l)a^ 



„i(2i/-l)a 



e "^^e^) 



[e^^\,^t))^'-^ Y: W{C) 

C|(a;,t)G7 

Again the same miracle happens in that ze^*"" — w. Note that the conservation law for ei is the obvious 
conservation of the flux of 7 through a plaquette. 

Finally, we define 



e['\x,t), (x,t) eZx (Z+i) 
\e"^er(x,t), (x,t) e (Z+ i) X Z, 



and we have 



h{x,t)) = ^ J2 W{C)e- 



-i9(C) 

z 

C|(a;,t)G7 

with the discrete holomorphicity equation 

(27) (j)i[x,t- 1/2) + e*("-"Vi(x + 1/2, i) - e'("-"Vi(a; - 1/2, t) - 0i(x,i + 1/2) = 0. 
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Note that since Ei commutes with each loop plaquette separately, cf. ([2T|) . the equation above is actually 
valid for each individual configuration of the loop model. 

Similarly, insertion of Ei leads to the antiholomorphic observable 

c\(x.t)ef 



In all cases of observables, we observe that the different Borel subalgebras {E operators versus E operators) 
correspond to different chiralities from the point of view of Conformal Field Theory. 



4.1.3. A remark on local observables. Note that the model also possesses a local observable, although its 
meaning is not transparent in the loop language. Write = q^'^ , so that in the representation we use, 
T^z{Hi) = —■Kz{Hq) = CT^. The coproduct of Hi is the usual Lie algebra coproduct A7?i = Hi ® 1 + 1 ® Hi, 



so that the corresponding current ft,'^^ = = —Iiq and /i*^*^ = — —Hq does not carry a "tail", i.e., 
is local. In the vertex language, (ft,(^^(a;, t)) simply gives the average orientation of the arrow sitting on the 
edge {x, t), and similarly for h'-^K This observable should not be confused with the flux observable associated 
to El or El. 



4.2. Application to the dilute loop model. 



4.2.1. Loop observables associated to Eq and Eq. Consider the dilute Temperley-Lieb loop model on the 
lattice described in § 12.11 with empty boundary conditions on all sides, except for one boundary site which 
is pointing inwards. The reason for this choice is that in the representation considered in § 12.51 the operator 
Eq turns a spin t into 0, or a into 4.- Therefore, in the loop model, inserting EQ{x,t) forces the path 7 to 
go from the boundary defect to {x,t). For definiteness, we place the boundary defect at the bottom, say at 



(xd, 1/2). Consider firstly the case where the operator sits on a horizontal edge. 




For the open path 7, we get the following factors. Denote the total winding angle of 7 by 6. It is an 
integer multiple of tt: 9 = kn (however we do not know anything about the parity of k this time because in 
the dilute model the w-tiles are parity-changing). Then e^ ^ gets a factor e*''^ from the local turns. Also, it is 
obvious that this loop crosses the tail (algebraically) [fc/2j times. This produces a factor of g^L^-'/^J from the 
tail (since the tail is comprised of To operators). Also, the terms with k odd get a factor q from the matrix 
element of Eq . The final expression for the time-component of the current eg*^ is 

/ \ 



ix,t) 



(p{q)w^~ 



E 

C|7:(a:<i,l/2)- 
. k even 



ix,t) 



E 

Ch:{xa.l/2)- 
k odd 



q e 



i(f-i)(fe-l)7r 



<x.t) 



W{C) 



DISCRETE HOLOMORPHICITY AND QUANTIZED AFFINE ALGEBRAS 



17 



which can be combined into a single summation 



Z 



Let us move on to the case where the operator sits on a vertical edge, Cq 



(^). 




We now have d — kn — a with /c e Z, and as in the previous case, the factor from the tail is g^L'^/zj ^ Putting 
everything together, we can again combine the two parities of k into a single summation 



C\'r:{xa,l/2)-^{x,t) 



So far the value of £ (the constant introduced in the evaluation representation tt^) has been left arbitrary. 

(x) 

In order to recover a discrete holomorphicity condition, we need the above formal expressions for Cg and 
Cq*'* to differ by a factor of e'" , and hence we require that 



(28) 

Recalling that (z/w) 



X e 



2£ 



one can check ([28]) is satisfied by choosing £ — 



We now define 



M^,t) :=^(g)-iu;^-i^°*^%/^' ^""'^^ ^ (Z,Z+i) 



C|7:(a;d,l/2)-^(2:,t) 

This observable satisfies the discrete holomorphicity condition: 

(29) (j)o{x, t ~ 1/2) + e'^^'-'^Uoix + 1/2, t) - e^'^^'-'^Uoix - 1/2, t) ~ M^, t + 1/2) . 

As in the dense case, choosing Eq instead of Eq would lead to the complex conjugate (f>o. 
Remark: If one defines ly ^ 5 ^ ^' i then we get 



Ch:ixa.l/2)-^{x.t) 



and X = z/w = e^^'^'^ which gives the same discretely holomorphic observable as in [T^, but a different 
relationship between the ratio of spectral parameters x and the angle a. 
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-id 



4.2.2. Loop observables associated to Ei and Ei. The operator Ei takes a spin 4, to t- Hence, as in the dense 
model, we should consider two boundary defects, keeping the rest of the boundary empty for example. The 
insertion of operator Ei at a given point then selects loop configurations such that the open path 7 that 
starts/ends at the boundary defects passes through this point. Let us now compute explicitly the resulting 
factors. Firstly there is a factor of w'^^ . For the winding of the 7 from the boundary to the point in the bulk 
(and from the point in the bulk to the boundary), we obtain a total factor of e"^*"^^, where 9 = kn. Finally, 
for the crossings of the tail (which is now composed of Ti matrices), we obtain a factor of g^'^ = g^i^"-^)^ , 
Putting everything together, we have 

2£ 

[e?{x,t))^^ Y: W{C)e- 

C|(a;,t)G7 

Notice that, in contrast to the dense case, here e~'^ is not independent of the configuration since the parity 
of k is not fixed. 

The calculation is completely analogous for ef\ Firstly there is a factor of z^^ . For the winding of 7 we 
obtain a total factor of e^^*"^, where 9 = kn — a. Finally, for the crossings of the tail we obtain a factor of 
q4k _ gi(2iy-i)(e+a) pitting everything together gives 

r)--^ — E mc)e-^ 

c|(x,t)e7 

Using = 22£g2jiyQ^ define a function on the edges of the lattice 

'e'i\x,t), {x,t) eZx {Z+ ^) 



e"'e["'\x,t), (x,t) G (Z+ i) X Z, 



0i(x,t) := w 

{Mx,t))^^ E WiC)e-^', 
c|(x,t)e7 

which satisfies the discrete holomorphicity condition 

Mx,t- 1/2) + e*("-"Vi(:E + 1/2, t) ~ e'("-"Vi(a; " 1/2,0 - Mx,t + 1/2) = 0. 
Again, using Ei we obtain the antiholomorphic counterpart (/>i. 

5. Vertex models with integrable boundaries 

5.1. Coideal subalgebras and i^T-matrices. Following the approach of [13], a left (resp. right) boundary 
quantized affinc algebra B can be defined as a subalgebra of U with the left (resp. right) coideal property 
A : i? — i> i? (x) [/ (resp. A : B ^ U ® B). If Vz is irreducible as a i? module, then the left (resp. right) 
boundary reflection matrix K£{z) : V^-i Vz (resp. Kr{z) : Vz — Vz-i) is the solution, unique up to an 
overall normalization, of the linear relation 

(30) Ki{z)TTz-i{Y)=7rz{Y)Kfiz) (resp. if,(z)^,(r) = 7r,-i (y)i^,(z)) 

for all Y ^ B. We take the following as the graphical representation of Kg{z) and Kr{z): 



Ki{z) = ( and Kr{z) 





where we recall that the arrows represent the flow of time as one reads algebraic expressions from right to 
left. 
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From this, Sklyanin [26] defined a "double- row transfer matrix" T2: 

T2(z; Wl, . . .,wl) = Tro {K,.^q{z)Rol{z/wl) ■ . ■ Roi{z/wi)Ki,o{z)Rio{zwi) . 
wliicli lias the graphical representation 



■ Rlo{zwl)) 



T2{z-,wi,...,wl) 



Wl 



By vertically stacking M double-row transfer matrices T2(zi; wi, . . . , wl), 1 < « < one builds a rectan- 
gular lattice of width L and height 2M, whose left/right boundary conditions are fixed a priori. Each T2 is 
an operator acting in V^^ ® ■ ■ ■ ® V^^ . 

5.2. Current conservation at the boundary. We now extend the formalism of [6] to the case of an 
integrable boundary. For simplicity, we treat only the case of the left boundary in full detail, but analogous 
relations can also be written for right boundaries. 

Suppose we are given a left coideal subalgebra C U, i.e., a subalgebra with the additional property 
A{Bi) (Z Bi® J/Q Assume also that we have a left boundary reflection matrix Ki{z) which satisfies for 
all Y G Bi. Provided that Ja G i?£, we have Kt{z)-Kz-i{Ja) = nz{Ja)Kg{z), which is represented graphically 

by 



(31) 



We have included tails of operators, but since with our conventions they extend to the left of the operator 
insertion, they never cross any lines and therefore play no role in equation (j31|) . This equation expresses the 
local conservation of the current at the left boundary. 

Similarly if 9a'' G Bi, we have Kg{z)'Kz-i{Qa'') = TTz{Qa^)Ki{z), which has the graphical representation 

Z n.u/Z 



(32) 




Analogous relations apply if Ja and 8a^ are elements of a right coideal subalgebra Br with right boundary 
reflection matrix Kr{z): 



(33) 



z^b 



and 



The result of the preceding equations is that we now obtain current and charge conservation laws which 
are exact, rather than correct up to boundary terms, as they were in 12.31 For example, current conservation 



In view of the simple coproduct relations l l3aH3bl l. one way to produce such a Bg is to have it generated by appropriate 
combinations of the elements J a and ©a''- We will discuss examples of such left coideal subalgebras in 15.51 
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([7|) is now automatic, since the assumption that the tail commutes with the left boundary is implicit in the 
condition Qa^ £ Bi. 

Similarly, by using equations ([5a|) and ([3T |) - ([33l) . we find that 





















< 

















































- ->- 






< 





































Graphically, this says that any charge built from JajQj' G Bi,Br commutes with the double- row transfer 
matrix; or equivalently, such a charge is conserved: 

(34) J,(i-l/2) = J,(t + 3/2). 

5.3. Light-cone lattice. Let us now consider a specialization of the parameters Wi on the vertical lines in 
our lattice. We assume L to be even and choose 

-1 



odd — even Z 



1<J<L. 



Assuming that R{1) oc P, which is the case for both U = UqiAl"') and U = Uq{A[^^), this causes every 
second i?-matrix present in the lattice to degenerate into a P-matrix. The resulting "light-cone" lattice has 
half the number of vertices, with the remaining ones being rotated by 45 degrees. The double-row transfer 
matrix becomes 




-2 — 



where we have absorbed an i?-matrix in the right boundary; denote the resulting boundary operator Kr{z) 
(we skip the details since we shall focus on the left boundary in what follows). All lines are now oriented 
upwards, so that we omit orientation arrows henceforth. Equivalently, T2 — T'^T", where 

i 

To :=Y[R2^+l{z^) = 

i 

Ri — Pi,i+ii?i,i+i is the i?-matrix acting on sites i, i -f 1 with an additional permutation P of factors of the 
tensor product, and Ki{z) (resp. K.r{z)) acts on the first (resp. last) factor of the tensor product. 

Since the light-cone lattice is obtained as a special case of the double-row, two-boundary lattice, all 
previous results continue to apply; they only need to be transcribed into the new orientation. The analogue 
of relation (I5ap is the commutation of R with the action of Ja, that is: 













(35) 



The 45 degree rotation makes "space" (resp. "time") lines go south-west to north-east (resp. south-east to 
north-west). The corresponding currents are simply 



Ja{x) = ^ 
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(or opposite tilt of the lines depending on parity of t), where x € Z and the former distinction between 
"time" and "space" components is uneccessary. The same conservation equation Q holds as a consequence 
of and the fact that 8q'' G Bi. Rewritten in the light-cone approach, it becomes: 

(36) ja{x,t)+ja{x+l,t)=ja{x,t+l)+ja{x+l,t+l), {x , t) e Z"^ , X + t = (mod 2) . 

Similarly, pip can be rewritten algebraically as 

(37) ja{l,t)^ja{l,t + l), t = (mod 2). 
There is a right boundary analogue: 

(38) jaiL,t)^ja{L,t + l), t = (mod 2). 
The charge J a is now the obvious sum: J a ~ 'l2x=ija{x), or graphically, 



(39) Ja - a 
and as a direct consequence of (|36H38l) . commutes with both To and To, i.e., 

(40) 3a{t)=3a{t+l). 

As already pointed out, we are mainly concerned with local current conservation rather than conservation 
of the charge; in particular the local conservation at the left boundary ([37]) only requires integrability at the 
left boundary. 

5.4. Adjoint action on the light-cone lattice. The adjoint action of an element Ja of U on another 
element Jt is defined as in ([9]). On the light-cone lattice, and in the case of an operator Ja commuting with 
ii'-matrices on all boundaries, the expression of can be greatly simplified. We illustrate this by a 

series of pictures in which, in preparation for the switch to loop models, we use dual lattice pictures, cf. (|16p 



in the bulk, and at the boundary, Ki = 



and similarly for the other boundaries. 



Using the conservation equation (j40|) . the contour defining ^abb*"*] can be deformed to follow the top and 
bottom boundaries, i.e., 
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We now suppose that top and bottom boundary conditions are also integrable, such that Ja commutes with 
the corresponding if-matrices, except at some boundary defects which sit say on the bottom boundary. Then 
the top part of the contour can be moved through the top row of iir-matrices, using again the boundary 
conservation relations. Indeed, the sum over the two points on each triangle along the top boundary gives a 
zero contribution!^ On the bottom part, the contour is reduced to an arch enclosing the defects: 



Translating this final picture into algebraic form, we obtain 

(41) {Aa[jb{xM = {[Ja{xd.Q)+%{xd + lMib{x.t)) , 

where we assume that there are two defects at adjacent locations Xd and + 1, and recalling that ja denotes 
a non-local current with a "right" tail, c/. 



5.5. Coideal subalgebras for quantized afRne algebras. Here we give examples of (left) coideal subal- 
gebras and boundary reflection matrices for the quantized affine algebras that interest us. Since we do not 
discuss analogous right boundary results, we omit the subscript "f" from all subsequent equations. 

5.5.1. The Uq{A^^^) boundary algebra. The choice of coideal subalgebra B wc shall consider in this paper is 
generated by 

{To ,Ti,Q:^Ei+rEo,Q Ei + rEo} , 



^To see this, we must pay attention to a crucial sign issue - our currents are associated to edges which are oriented upwards, 
so current conservation at the top/bottom boundaries must be accompanied by a sign in one of the terms. It is precisely this 
sign which causes the pairwise cancellation for each triangle. 
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where r is a real parameter. The left coideal property is satisfied because 

A{Q) ^Q(g,l + Ti(^Ei+To®rEo and A{Q) ^ Q I + Ti (g, Ei + Tq rEo 
are both elements oi B ^U. After a choice of normalization the solution of ([5(1)) is [H] 



(42) 



K{z) = 



z + rz ^ 







z ^ + rz 



Note that K{z) is diagonal as a consequence of the fact that Tq and Ti are elements of B. 

(2) 

5.5.2. The Uq{A2 ) boundary algebra. The boundary algebra B in this case is taken to be that generated 
by {To,Ti,Ei,Ei,Q,Q}, where 

Q := [EuEo]q-^+rEo, Q := [Eu Eo]g^ ^ rq^ Eq . 

in which we use the notation [a, foj^ = ah — xba. 

The left coideal property A(i3) E B ^ U is easy to check. With a choice of normalization, the solution 
of ([30]) in the case when r is fixed to be r = ±iq^^ reads [H [22] 

z'^'^iz-^ +rz) 

(43) Kiz) = ■ 



z + rz 






z^^^{z^^ + rz) 
(1)^ 



For definiteness, we shall henceforth take the root r ~ iq ^ . As in the Uq{Ai ') model case , K{z) is diagonal 

because Tq and Ti are elements of B. In contrast to the Uq{A\^^) case, there is no solution of (pO)) for general 
values of the parameter r. 



6. InTEGRABLE BOUNDARIES FOR LOOP MODELS 

In this section we repeat the ideas of §[3] to introduce boundary tiles into the dense and dilute Temperley- 
Licb models. In complete analogy with § [31 the corresponding K matrices (|^ and P5)) are recovered as 
linear combinations of the boundary tiles. We use the light-cone approach of § 15. 3[ with an angle of a on 
the lattice, i.e., 



R 




K = 



and similarly for other boundaries. 

6.1. The boundary dense Temperley Lieb model and the J7g(Aj^'') vertex model. Let us define a 
boundary loop model by introducing an additional weight 1 boundary plaquette, as follows: 
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and by assigning weight r^") — — (e*'^(^'^^"4) + g-'i'(27r-n?)') \oo\i that that passes n times through a 

boundary. Thus we can view the boundary as introducing a deficit angle of ^. 

Again, we can turn this boundary weight into that of a vertex model by viewing it as an operator from 
— )■ in a S-N direction, resulting in a boundary weight 

e"'("-«) ) ■ 

Since the spectral parameter z is related to the angle a by z = e"*''", this boundary matrix coincides with 
Kiyz) of dlilQ if we take 

g2»^« ^ (l + rz-2)/(i + rz2). 

Clearly the r = case corresponds to a zero deficit angle ^ = in which case the boundary TL plaquette 
becomes the single plaquette with weight one which we call "free boundary conditions" and denote by 




f 2") 

6.2. The boundary dilute Temperley Lieb model and the ) vertex model. We introduce 

two additional boundary plaquette configurations with associated weights p and k, as follows: 



( 



Interpreting the plaquettes as operators as above yields a boundary reflection matrix 
(45) K = 



which is equal to the matrix if (z) of equation (|43)) if p = z + ^ and k ~ z ^ + ^z, using also the 

,2£ 



— %VCL 











p 














7. Non-local currents and boundary discrete holomorphicity 

7.1. Application to the dense loop model. For convenience, in this section we shall use exclusively the 
light-cone orientation of the lattice. We consider loop configurations on the lattice which contain a single 
open path 7. For simplicity, we assume that the ends of 7 are situated next to each other, as follows: 



^In fact the two X-matrices coincide after renormalizing I I44II by ((1 + rz ■^)(1 -|- rz^))^/^. Since this only produces a global 
factor in the observables we will consider, we omit this normalization for simplicity. 
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As we did in the case of trivial boundary conditions, we consider observables which are constructed by 
requiring that the open loop goes through a certain point {x,t), either in the bulk or on the boundary of 
the lattice. Because the i^-matrices are diagonal the correct way to obtain such observables is, as before, 
to insert a local operator (say Eq, complete with its tail) at since all lattice configurations vanish for 

which Eq is situated on a closed loop: 




En 



= 0. 



7.1.1. Loop observables in the bulk. Let us repeat the analysis of the observables in the bulk, but now in 
the presence of non-trivial boundaries and on the light-cone lattice. Insert eg at the point {x + l,t) £ Z^, 
X + t = (mod 2), in the bulk of the lattice. As mentioned above, the only configurations which survive are 
those for which the open loop goes through {x + l,t). The contribution of the open path to the weight can 
be found in a similar way as before as g^'^i^^-'^+^i) qf^ ^ where 6 = kn is the angle formed by the left-incoming 
portion of the open loop (if we treat the boundary tiles on equal footing with bulk tiles) , and n the number 
of contacts of the left portion of this loop with the boundary minus that of the right portion. So we find 



{eQ{x + l,t)) 



Z 



■ E 

C\{x+l,t)£', 



W{C) e''^'^"-^^^ e^'"^ , 



Similarly, repeating for eQ{x,t) with a; + 1 = (mod 2), we obtain 



(eo(x,i)) 



z e 



C|(2:,t)67 



The local conservation law (|36l) with ja — eg is 
(46) eo{x,t)+eQ{x + l,t) = eo{x,t + I) + eo{x + l,t + I) , 



{x,t) eZ^ x + t = (mod 2). 



This holds in the bulk because the tail, which is comprised of Tq operators, commutes with the left X-matrix 
(since To G B, as explained in § I5.5.ip . Therefore, using analogous arguments to those of § 14. 1[ we define 
the function 



4>o{x, t) :— ze*' 



eo{x,t), x + t = l (mod 2) 
e'°'eoix,t), x + t = (mod 2) 
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and in view of the fact that z = e we have 

(47) (0o(^,i))^i W^(C) e^(4''-i)^e™''« . 

c|(x,t)e7 

Applying pS)) to this observable, we find that it satisfies 

e*('^^"Vo(2;, t) + (t>o{x + f , t) - e'^''-°'Uo{x + 1,^ + 1)- M^, t + 1) ^ , 
which is discrete holomorphicity on the light-cone lattice. 

The observable corresponding to insertion of Ei gets modified in an analogous way; namely 

\ei{x,t), x + t = l (mod 2) 



(/>i(a;, t) := ze 



e'°'ei{x,t), x + t = (mod 2) 



gives rise to the observable 

(48) (0i(x,t)) = l J2 M^(C) e-'*' e-™^« . 

This is the flux observable in the presence of a non-trivial boundary and it satisfies the discrete holomorphicity 
equation 

e'("~"Vi {x, t) + ^i{x + 1, t) - e'^'^'^Uiix + l,t + l)-<j>i(x,t + l) = 0. 

7.1.2. Loop observables at the boundary. By boundary we mean in what follows the left boundary. Since 
r 7^ 0, the operators Eo,Ei,Eo,Ei are not elements of the coideal B, and hence they are not conserved 
separately at the boundary (with trivial boundaries, Ei and Ei were conserved), and in particular their 
associated charge will not be conserved. 

We now consider the combinations Q — Ei + rEo and Q = Ei + tEq, which are in B. Using ([57]) in the 
case ja = ei + reo we find that 

ei(l,i)-|-reo(l,t) = ei(l,t + l)+reo(l,i-Kl), t^O (mod 2) , 

which can be translated into the following equation for the observables (f>i and (f>o: 

.2- Vi(l, t) + rzM^, t) = e-*"z-i0i(l, t + l) + e'^rzM^, t + 1) . 

This is neither a holomorphicity nor an antiholomorphicity condition because we are mixing operators from 
the two chiralities. However, by taking the real part (or equivalently, summing this identity and the one 
satisfied by the conjugate observable Q), one finds that ip := z~^{4)i -\- r(j)o) satisfies 



Re 



V'(l,t) + e'('^"")V^(l,t + l) 



0, 



which is a boundary discrete holomorphicity condition around the plaquette 

C a t = (mod 2). 
,•(1,*) 

Remark: At the left boundary the tails (which are on the left) disappear and therefore the two observables 
01 and 00 are the same up to a constant. This can be seen more explicitly in the fact that there cannot be 
any winding at the boundary, so the angle 9 in ()47|48p is fixed and independent of the configuration. 
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7.2. Application to the dilute loop model. We consider the dilute loop model on the light-cone lattice, 
with as in the dense case possible defects located next to each other on the bottom row. All the observables 
that we consider force a line to go from the insertion point to the boundary defect, so that a typical 
configuration is as depicted below. 



7.2.1. Loop observables in the bulk. In contrast with the dense case, the X-matrix (j43[) does not introduce any 
orientation-dependent phase factor, but simply a relative weight for contacts with the boundary. Therefore, 
the analysis of observables Eq and Ei in the bulk is unchanged, and we do not repeat it here. 

A natural question is whether one can use the adjoint action to build new currents. We consider the 
element P G UqiA^^^) given by P = a,dEi {Eq). Explicitly, 

P — EiEq — TiEqT^ ^Ei — EiEq - q^^EqEi , 



and hence we can use the treatment of § 15.41 to express the corresponding lattice observable p. From ()4ip , 
we have 

{p{x,t)) = {Ae, [eo{x,t)]) = {[ei{xd,0) + ei{xd + l,0)]eo{x,t)) . 

Since Ei can only flip a state I to a. state t, the non-zero terms correspond to a pair of defects (J,, 0) or (0, 1). 
Summing over these two possibilities, we get: {p{x,t)) = — g~^(z^^ -I- {eo{x,t)). Hence, the operator P 
does not lead to a new holomorphic observable; it simply produces (f>o, up to a multiplicative constant. 

Equivalently, let us define 

^(x t) ■= z^~^ V^^'^^' a;-ft = l (mod 2) (i fx ■= z^" V x + t^l (mod 2) 

^ \e'°'p{x,t), x + t^O (mod 2), ^o(x, ) . z a; + 1 (mod 2) . 

Then 

(49) {i{x, t)) = c (0o(a;, t)) , where c = -q^^{z^^ + z^^O- 

Since it is only possible to relate f and 0o a-s expectation values (•••), the constant c is dependent on the 
choice of boundary conditions; it would differ, had we chosen alternative boundaries to the ones shown in 
the figure above. 

7.2.2. Loop observables at the boundary. We are now in a position to construct an observable which satisfies 
discrete holomorphicity at the left boundary. We consider the operator Q = P + iq~-^Eo (with P as in the 
previous section), which commutes with the i^-matrix in the sense of ([30|) . cf. §5.5.21 
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From p7[) with ja = p + iq ^sq we obtain 

p(l, t) + iq~^eoil,t) = p{l,t + 1) + iq-^eo{l, t + 1)., 
which can be translated into an equation for the observables ^ and 



zi-*^(l,t) + zg-iz^-Vo(l,t) = e- 
Taking the real part of this equation, we obtain 



'z^-'^il, t + l) + e^^iq-h'-^Mh t + 1) 



Re 



1/^(1, <) + e'("-"V(l,i+l) 



= 0, 



where -0 = 2^ ^(^ ~ iq'Po)- Now taking into account (|49p . we find that in the equation above one can replace 
ip with e*^0o (A being some phase, dependent on our choice of boundary conditions, which we do not write 
explicitly) . 



8. The continuum limit 



8.1. Dense loops. In this section, we identify the operators corresponding to the lattice holomorphic ob- 
servables, as holomorphic currents in the CFT describing the continuum limit. 

8.1.1. Scaling theory. Let us first briefly review the Coulomb gas construction [23l|T2]. We define a height 
function ^{x + 1/2, t + 1/2) on the dual lattice, such that the oriented loops are the contour lines of $, and 
with the convention that the values of 4> across a contour line differ by n. In the continuum limit, the height 
function is then subject to a Gaussian distribution: 

(50) S[^] = Y j i^^fdxdt , where g^l-2v. 

Note that the coupling constant spans the interval < p < 1. In what follows, the model is defined on 
a cylinder of even circumference L, and the axis of the cylinder is in the time direction. Because of the 
definition of $ by local increments, $ can be discontinuous along the circumference: 

^(i + .T + l/2,t + 1/2) - + l/2,t + 1/2) = 27rm , m e Z. 

Thus, the height function should be considered as living on a circle: 

$ = $ + 27r . 

Also, the local Boltzmann weights associated to loop turns ensure that the closed loops get the correct 
weight r = 2cos(27rt/), except for the non-contractible loops: these loops have a vanishing total winding, 
and thus get a weight t = 2. To restore the correct weight, one introduces a seam in the time direction, 
such that every right (resp. left) arrow crossing the seam gets a weight e"" (resp. e"*'^"). The weight of 
non-contractible loops becomes r = 2cos7ra, and one sets a := 21/ to get t = t. 

8.1.2. Operator content. To recover the full-plane geometry, we use the complex coordinates 

In this setting, the seam described above goes from the origin to infinity, and amounts to introducing a pair of 
vertex operators e'"*(°°)e~'"*^°^. More generally, if we decompose the height field as z) = (p{z) + (p{z), 
the primary operators are of the form = ^^ip-v+P-v) ^ which we write as 

The first factor is only well-defined if n := {fi + /U,)/2 is an integer, which is called the electric charge. The 
average value of the second factor is e**"^', where $ci = — M)(log^ — log2)/(4g). The discontinuity of 
$ci around the origin is (5$ci — 2t: x {fi — jl)/{2g), and hence the number m :— {fi — Jl)/{2g) must be an 
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integer, which is called the magnetic charge. Since the conformal weight of the chiral vertex operator e*^*^ 
is hf^ = fi{pL + 2a)/{4:g), we get for O^,,^,: 

^ {n + a + mgY — o? ^ (n + a — mgY — o? 

with n^m £ 1? . In this context, a. appears as a background electric charge. Part of the above spectrum fits 
in the Kac table for minimal models: 



4<? 

with integer r, s. 



8.1.3. Scaling limit of the lattice observables. We wih now show that the discrete (anti-)holomorphic observ- 
ables discussed in this paper scale to operators of the form 0^,p. Let us first consider the two-point function 
associated to Eq: 

C|7i,72:tu-*-^ 

where the sum is on loop configurations with two open oriented paths 71, 72 going from w to z, and Q is the 
winding angle of each path. In terms of the height model, this correlation function includes magnetic defects 
of charges —1 and +1 at z and if, respectively. Moreover, the winding angle is given \iy d — $(2:, z) — ^(w, w), 
and hence the factor e'^**""^)^ corresponds to an electric operator of charge n ~f a = 4;^ — 1 at z, and opposite 
charge at w. These charges yield the values \x — —2g and /l = 0, and thus we identify: 

with conformal dimensions h = 2g — 1 = /113 and h = 0. By reversing the arrows, we obtain 

00 (z) = e-^^f-^t^") . 



Similarly, the two-point function associated to Ei is: 

C|7i,72:'"'— >2 

which has charges m = 1 and n -f a = 1 , and hence we get 

01 (z) = e^ '^g^i-) ^ 01 (^) ^ ^-i^g^i^) . 

The holomorphic current 0i has conformal dimensions h = 1 and h — Q: this is the "screening operator". 

We now turn to the lattice operator associated to the diagonal generators Hi. Since Hi oc u^, it simply 
measures the local orientation of loops. The increment of $ across an up (resp. down) arrow is tt (resp. — tt), 
and thus one has adx^ — tt/i*^*-', and likewise adt^ — —'Kh^^\ where a is the lattice mesh size. Using the 
complex coordinates w = t + ix and w — t — ix, we identify the chiral currents: 

The conservation law found in § 14.1.31 corresponds in the continuum limit to the conservation of non-chiral 
current d*^ — d*{tp + (f), or in components, e'^'^di,^, which ensures local well-definedness of $. 

The above results can be summarised in the following diagram: 
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d 




In this figure, the horizontal axis is the U{1) charge a^, and the vertical axis is the gradation d in the 
evaluation representation of A'{\ We note the similarity with the discussion of nonlocal charges in the 
(ultra-violet limit of the) sine-Gordon model in [7]. A notable difference is the choice of gradation, which 
has different origins in the two situations. 

8.2. Dilute loops. The mapping to a compactified free boson CFT (|50p also holds in the dilute case, up to 
minor adaptations. We keep the convention that the height function $ has jumps of ±7r across an oriented 
loop. Since empty edges are now allowed, this means that the discontinuities of $ along the circumference 
of the cylinder are now multiples of tt, and one should set $ = $ + tt. So we can keep the same notations 
as in the previous section, except that the allowed electromagnetic charges become: m G Z/2 and n G 2Z. 
Moreover, ly is now chosen in the interval [—1/2, 0], and we have 1 < i? < 2. 

The two-point function for 0o is: 

where the sum is on loop configurations with one open oriented path 7 going from w to z, and 9 is the 
winding angle of 7. Since this is a one-leg defect and 6 = 2[^{z,z) — ^{Wjw)], the corresponding charges 
are m = —1/2 and n + a = ?>v — 1/2. The "flux observables" </>! and 0i are the same as in the dense model. 
Thus we obtain 

(j,^ (z) = e^'svi^-^ , 01 (z) = e^'^^^'^ , 

and the conformal dimension for (j)o is h — {3g — 2)/4 = hi2, whereas for it is ft, = 1. Finally, the diagonal 
operators /i^^'*) relate to dw(p and dw(f. 

9. Conclusions 

In this paper, we have described a general procedure to obtain discretely holomorphic observables out 
of nonlocal currents in quantum integrable lattice models. We have shown in several examples how these 
observables are naturally expressed in terms of loop models. We have identified them in the continuum limit, 
connecting to Conformal Field Theory. 
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It should be noted that in CFT the conserved currents always come in pairs: a current and its dual 
current — e^^j". Only the two conservation laws combined imply separation of chiralities, and therefore 
existence of holoniorphic observables. Here we only base our analysis on a single conservation law for each 
observable, hence a "weak" discrete holomorphicity condition - the dual equation is missing. This absence 
can be traced to the step in which we identify the two components (say, time and space) of our current 
as a single function corresponding to the observable. This step would require additional justification in 
order to proceed with a rigorous proof of the conformal limit|j The fact that in all cases, our would-be 
chiral observables have, in the loop language, a unifying definition on both vertical and horizontal edges is 
certainly a strong indication that such an identification is correct. 

This work opens the way to further study and interpretation of discrete holomorphic observables, in 
particular in the case of more general boundary conditions (as recently studied in Also, the application 

of this approach to off-critical models (see the treatment of the Ising model in [3S]) needs to be developed. 
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